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In this paper, we calculate the next to the leading order Casimir energy for real massive and 
massless scalar fields within \<p 4 theory, confined between two parallel plates with the Dirichlet 
boundary condition in two spatial dimensions. Our results are finite in both cases, in sharp contrast 
to the infinite result reported previously for the massless case. In this paper we use a renormalization 
procedure introduced earlier, which naturally incorporates the boundary conditions. As a result 
our radiative correction term is different from the previously calculated value. We further use a 
rcgularization procedure which help us to obtain the finite results without resorting to any analytic 
continuation techniques. 

PACS numbers: 



I. INTRODUCTION 



The Casimir effect can be observed in all systems with 
nontrivial boundary conditions (BCs) or background 
fields (e.g. solitons). During the last sixty years this ef- 
fect has been an important topic of research with ap- 
plications in many branches of physics [3, 0, S| . The 
static Casimir effect, first calculated in 1948 0], predicts 
an attraction between two perfectly conducting paral- 
lel plates, due to distortion of the electromagnetic vac- 
uum state (for a general review on the Casimir effect, see 
Refs. @, [a, |(|). Ten years later the first attempt to ob- 
serve this phenomena was made by Sparnaay |7j. Since 
then, many experimental investigations have measured 
precisely the Casimir force in different cases, such as two 
parallel plates [|, or a sphere in front of a plane [£$ . The 
majority of the theoretical investigations related to the 
zero order Casimir effects are for various fields, geome- 
tries, BCs [1, H3, EH EH, an d various dimensions. Some 
of the major approaches used are: the mode summa- 
tion method with a combination of the zeta function reg- 
ularization technique [H, E3 , Green's function formal- 
ism [iBK multiple-scattering expansions [l6j . heat-kernel 
series [171 ] . On the other hand, there exist many works on 
the first order and also second order radiative corrections 
to the Casimir energy for various cases [18|, [l9|, 120, |2l| . 
Some of the major approaches used for the radiative cor- 
rections to the Casimir effect are the phase shift of the 
scattering states [22j], or the replacement of the BCs by 
an appropriate potential term [23|, [24[ ■ The Casimir ef- 
fects have found many applications in physics. For ex- 
ample, the Casimir effect is the major contribution to 
the radiative correction to the mass of the solitons, and 
these corrections have been investigated in many pa- 
pers [H, US H3 . 
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The value of the Casimir energy has a complicated 
behavior as a function of the number of spatial dimen- 
sions, the type of fields, type of topology, and geometry. 
The case of even spatial dimensions is usually more com- 
plicated [12, HH, HH . Since many interesting condensed 
matter systems are well- approximated by two dimen- 
sional models, extracting finite results from the compli- 
cated divergencies, which usually plague such systems, 
is very important. The Casimir energies for scalar fields 
in even dimensions have been discussed for two parallel 
plates [2l[ , spheres [H, , and cylinders [TJ] . Some of 
those cases give divergent results, and some authors pre- 
scribe methods to extract a finite answer from those ex- 
pressions [l,E2,Eil]- However, those prescriptions are not 
universally applicable to all even dimensions, ft seems 
that even for the simplest case of a scalar field in 2 + 1 
dimensions, it is not clear how the divergences can be 
removed [a] ■ 

In this paper we calculate the first order radiative cor- 
rection to the Dirichlet Casimir energy for two infinite 
parallel plates for massive and massless scalar fields in 
two spatial dimensions. The problems mentioned before, 
give us extra motivation to utilize an alternative renor- 
malization program and regularization procedure for this 
problem. We have used these procedures to calculate this 
quantity in 1 + 1 and 3+1 dimensions [29|, ya]. As we shall 
see our procedure yields finite results for both massive 
and massless scalar fields which is different from the pre- 
viously reported one [2~l| . As a matter of fact the previ- 
ously reported result for the massless case is infinite [2l[ ■ 
It is worth mentioned that our finite result is obtained 
without any use of analytic continuation techniques due 
to our regularization procedure. The difference between 
our results and the previously reported one can be at- 
tributed to our alternative renormalization program. 

In our paper we combine two independent programs in 
order to calculate the radiative correction to the Casimir 
energy. First, we use an approach to the renormaliza- 
tion program which we believe to be systematic. The 
procedure to deduce the counterterms from the n-point 
functions in the renormalized perturbation theory is stan- 
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FIG. 1: The geometry of the two different configurations 
whose energies are to be compared. The labels al, 61, etc. de- 
note the appropriate sections in each configuration separated 
by plates. The left configuration is denoted by 'A' configura- 
tion, and the right one by 'B' configuration. 



dard and has been available for over half a century [3l| . 
We believe that all of the information about the nontriv- 
ial BCs or position dependent background fields should 
be carried by full set of the n-point functions. There- 
fore, all of the counterterms deduced from these n-point 
functions should also contain these information. Using 
this procedure we deduce the position dependent coun- 
terterms in our problem. We should mention that most 
of the authors use the free counterterms, by which we 
mean the ones that are relevant to the free cases with 
no nontrivial BCs, and are obviously position indepen- 
dent. However, the dependency of the counterterms on 
the distance between t he p lates has been noted in some 
references such as [HI, l33j . However, these authors use 
free counterterms in the space between the plates and 
place additional surface counterterms at the boundaries. 

Another important part of our calculation is using a 
method to remove the divergences without resorting to 
any analytic continuation. In fact, we subtract two dif- 
ferent configurations with similar nature. This subtrac- 
tion scheme is based on the Boyer's subtraction scheme 
and it can regularize the infinities and help us to remove 
them without using any analytic continuation I3JL This 
method has been used in many previous works |3l. l35l. |36| . 
We place the two infinite parallel plates (with distance 
a) within two other plates (with distance L > a). We 
then construct a similar configuration of plates with dis- 
tances L > b. We then subtract the Casimir energies of 
these two configurations. Finally in order to obtain the 
Casimir energy for the original configuration we let L 
and then 6 go to infinity. Therefore, the Casimir energy 
is now defined by, 
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where 
E A 



E a i + 2E a 2, 



Eb — Eti + 2Eb2 , (2) 



and E a i, E a %, E^i and E^ are the zero point energies of 
each region shown in Fig. ((T|). 



We have already used this subtraction scheme to calcu- 
late the leading order part of the Casimir energy for a real 
massive scalar field and its massless limit with Dirichlct 
BC for two infinite parallel plates in arbitrary dimensions 
in Refs. [29|, [3(|. Therefore, in this paper we only report 
its final result for two spatial dimensions. We obtain, 
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This expression for the leading order of the Casimir 
energy of a massive scalar field with Dirichlet BC in 
two spatial dimensions is the same as that reported in 
Refs. 3, jlfJl . However, contrary to the methods used in 
Refs. 3l- 1 1 CM] . this expression is obtained without using any 
analytic continuation techniques. Two important limits 
should be considered at this stage. First is the small mass 
limit, m — * 0, and Eq. (J3]) becomes, 
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where £(s) denotes the zeta function. This expression 
for the leading order Casimir energy of a massless scalar 
field is also the same as reported in Refs. [!, [l(| ■ Second 
is the large mass limit, ma 1, and Eq. ([3]) becomes, 
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In fact, in this limit the value of the Casimir energy 
decreases exponentially with increasing ma, and this is 
again the same as the previously reported result [!, • 
In section 2, we first calculate the first order radiative 
correction to the Casimir energy for this problem. We 
then plot all of the result for the massive and massless 
cases. In section 3, we summarize and discuss our results. 



II. FIRST-ORDER RADIATIVE CORRECTION 

In this section we first calculate the leading order ra- 
diative correction to the Casimir energy for a massive 
scalar field within X<p 4 theory with Dirichlet BC in 2 + 1 
dimensions using the renormalized perturbation theory. 
As mentioned in the Introduction and also Refs. [29 . |3CJ | , 
the counterterms are computed from the appropriate n- 
point functions which, in the presence of the nontrivial 
BCs, are naturally position dependent. The renormaliza- 
tion procedure, the deduction of the counterterms, and 
the final general form of the first order correction to the 
Casimir energy for each region have been completely dis- 
cussed in Refs. [lij [3(j. Therefore, in this paper we use 
only the conclusions: The general expression for the first 
order radiative correction term to the Casimir energy is, 
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where G a i(x,x') is the propagator of a real scalar field 
in region al in two spatial dimensions. After the usual 
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wick rotation, the expression for the Green's function or space is, 
the propagator G a \ (x, x') in three-dimensional Euclidean 
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All of the integrals in Eq. ([5]) are logarithmically diver- 
gent, and we make them dimensionless by multiplying 
appropriate factors of a. Then, we use cutoff regulariza- 
tion for each integral, and expand the results in the limit 
in which the cutoffs go to infinity as follows, 
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Using Eq. (J9j) for each integral in Eq. (JSj> we obtain, 
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where uj' 2 x n = (nir) 2 + m 2 a 2 , and A al is a cutoff in the upper limit of the integrals in Eq. ©. The terms related to 
other regions in Fig. JT]) are calculated analogously. Now, for the calculation of the Casimir energy in Eq. ([I]), we have 
four similar terms which should be subtracted from each other. By appropriately adjusting each cutoff A a i, A a 2, Af,i 
and A/,2, all of the infinities cancel due to our box subtraction scheme. We obtain, 
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Now we can use the Abel-Plana Summation Formula (APSF) which basically reduces the summations into integrations 
as follows, 
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If the summation starts from n = 0, the sign of the first term becomes positive. Now by applying the APSF to all of 
the summations in Eq. (|11[) we obtain, 
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where 

U(x) = -(^-lnm 2 x 2 + [ dnln (n 2 ir 2 + m 2 x 2 ) + B 1 (x) ]+ ^- ( In 2 m 2 x 2 + j dn In 2 (n 2 Tr 2 + m 2 x 2 ) + B 2 {x) ] , 



J dn In (nV + m 2 x 2 ) + B Y (x) j+ 
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and B\(x) and B2{x) are the branch-cut terms of the 
Abel-Plana summation formula and are calculated in Ap- 
pendixlA"! Both of these two types of branch-cut terms 
are finite for m ^ 0. However, other integral terms 
which appear in Eq. Q13[) are divergent. At this stage our 
main purpose is to regularize these terms and show how 
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they cancel each other, again due to our box subtraction 
scheme. 

To regularize the integrals in Eq. (fT5|) . we set separate 
cutoffs, denoted again by As, for the upper limits of each 
integral. After the integrations, we expand the results in 
the limit A — > oo. Now, by appropriate adjustment of the 

I As, all of the divergent terms which depend on the cutoffs 

As, cancel in Eq. (| 13[) . due to our box subtraction scheme. 
Below, we present the details of these cancelations for 
both types of integrals. For the first type we have, 
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where in the first line we have used the following change of variable N = nit /ma. Therefore, only the finite terms 
{am, (L — a)m, bm, {L — b)m} remain for the first type of integrals. For the second type of integrals we have, 
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where /C(A) in the limit A — > oo is, 
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The first term in the right hand side of Eq. (|16p in line 
three is divergent. The second integration is similar to 
the first type of integral terms which were considered 
above. The divergent terms in the sum of these two terms 
are removed by choosing appropriate adjustment of As 
and using the subtraction scheme indicated in Eq. (fT3)) . 
and only finite terms remain. The third term in the right 



hand side of Eq. (fTI)]) and its counterparts in the other 
regions are also divergent and their calculations are very 
difficult. However, if we let first the cutoffs go to infinity, 
one can show that they exactly cancel each other in the 
box subtraction scheme. Therefore, the only contribu- 
tions coming from this term is, 
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Using Eqs. (|13I15I18)) . we have 
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There many internal cancelations in the above expres- 
sions. After these cancelations only the branch-cut terms 
remain. By using the values of the branch-cut terms ob- 
tained in the Appendix^ we can write an explicit ex- 
pression for the lowest order radiative correction to the 
Casimir energy in terms of parameters m, a, ^^,b and 
^Y 1 - As stated in the Eq. first the limit L/b — > 00 
should be calculated and then b/a — > 00. In these limits 
all of the terms which depend on L and b disappear from 
our expression and only the terms which depend on the 
distance of the original plates (a) remain. Our final result 
is, 
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where 7 is the Euler-Mascheroni constant, and T(a, x) is 
the incomplete Gamma function. Our result differs from 
the previously reported result [2]| , since they use the free 
counterterms, and we have used the ones dictated by the 
Green's function appropriate for this problem. 

To calculate the massless limit, we go back to the orig- 
inal expression given in Eq. Q19p . The direct calculation 
of the massless case is extremely difficult. We use m as 
a regulator for this limit. However, multitude of difficul- 
ties appear. These difficulties are partly due to the fact 
that the branch-cut terms are also divergent in the limit 
m — > 0. Fortunately, there is no essential singularity and 
we obtain 
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FIG. 2: The first order radiative corrections to the Casimir 
energy for massive and massless scalar fields in two spatial 
dimensions are plotted as a function of the distance between 
the lines (a), within the \(f> 4 theory for A = 0.1. The numerical 
values for the plots have been multiplied by a factor of 100, 
in order to make their absolute values comparable to the zero 
order terms shown in Fig. @. In this figure we have shown the 
sequence of plots for m = {1, 0.1, 0.01, 0.001, 0}. It is apparent 
that the sequence of the massive cases converges rapidly to the 
massless case and there is an insignificant difference between 
the figures of the massive cases for m < 0.01, and the massless 
case. 



where in the second line we have used the small mass 
limit of B±, 



-2ma\ ™— >9 



81(a) = ln(l - e - zma ) 



ln(2ma) - ma + 0(m 2 ), (22) 



and used a suitable change of variables for B2 which 
leads to some cancelations, and we have ignored terms 
of 0(m 2 ). In the above expression all of the infinities 
cancel and we finally obtain the following finite result, 
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As shown in Fig. ([2|), the sequence of plots of the massive 
cases converges rapidly to the massless limit. It is obvi- 
ous that the massless limit is finite, exactly as we have 
obtained, and as expected on physical grounds. 

In Fig. all of the values for the zero order and the 
first order radiative correction to the Casimir energy for 
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FIG. 3: The leading term for the Casimir energy and its first 
order radiative correction (multiplied by a factor of 100) in 
two spatial dimensions, are plotted as a function of the dis- 
tance between the lines (a) for a massive (m = 1) and a mass- 
less scalar fields for A = 0.1 . The correction terms are always 
positive. 



parallel plates in two spatial dimensions has been cal- 
culated by a systematic approach to the renormalization 
program. This program automatically yields position de- 
pendent counterterms. Moreover, we used the Boyer's 
subtraction scheme which eliminates the need to use any 
analytic continuation techniques. The final results for the 
radiative correction of the Casimir energy for a massive 
and massless scalar fields are different from the reported 
results in the previous papers [2l| ■ We believe that this 
difference is due to the use of different renormalization 
programs. It is important to note that our result for the 
massless case is finite, in sharp contrast to the previously 
reported result [2l[ . 



a massive (m = I) and massless scalar fields are plotted. 
We should mention that the correction terms are posi- 
tive and their values arc approximately 100 times smaller 
than their zero-order counterparts. 

III. CONCLUSION 

In this paper, the first order radiative correction to 
the Casimir energy with Dirichlet BC for two infinite 



APPENDIX A: CALCULATION OF THE 
BRANCH-CUT TERMS 



In this Appendix we present the calculation of two 
types of branch-cut terms which appear in the calcula- 
tion of the first order radiative correction to the Casimir 
energy. We start with the first type of the branch-cut 
term which is denoted by B\ in the main text. We have 
for -Bi(a), 
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The values of other branch-cut terms can be easily written only by the appropriate replacement in the argument of 
B\. Analogous process is repeated for calculation of the second type of the branch-cut term. So, for i?2(a) we have, 
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This integral can not be performed in closed form. Ex- 
panding the denominator of the integrand we obtain, 
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where T(a,x) is the incomplete gamma function and in 
our case we have, 
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